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We establish by exact, nonperturbative methods a uni- 
versality for the correlation functions in Kraichnan's "rapid- 
change" model of a passively advected scalar field. We show 
that the solutions for separated points in the convective range 
of scales are unique and independent of the particular mech- 
anism of the scalar dissipation. Any non-universal depen- 
dences therefore must arise from the large length-scale fea- 
tures. The main step in the proof is to show that solutions of 
the model equations are unique even in the idealized case of 
zero diffusivity, under a very modest regularity requirement 
(square-integrability) . Within this regularity class the only 
zero-modes of the global many-body operators are shown to 
be trivial ones (i.e. constants). In a bounded domain of size 
L, with physical boundary conditions, the "ground-state en- 
ergy" is strictly positive and scales as L -7 with an exponent 
7 > 0. 

PACS numbers: 47.10.+g,47.27.-i,05.40.+j 



The Kolmogorov 1941 theory (K41) of fully devel- 
oped turbulence hypothesized the existence of universal 
statistics at so-called inertial-range scales L 3> £ 3> rj, 
where L is the integral length characteristic of the peak 
energy and 77 is the (Kolmogorov) microscale character- 
istic of the peak dissipation. According to the first and 
second similarity hypotheses of Kolmogorov jjj , both the 
limits L — > 00 and r\ — > should exist for the distribution 
function (PDF) of appropriate inertial-range variables, 
such as the velocity-differences 6ve(x) — v(x + t) — v(x), 
and depend only on the length-scale I and the mean dissi- 
pation e per mass. However, much evidence has accumu- 
lated subsequently that there is a nontrivial dependence 
of such PDFs on the ratio L/t This is associated to 
an increase in the fluctuations, or "intermittency," of the 
inertial-range variable for increasing ratio L/£, so that 
the first limit L —> 00 appears to lead to diverging mo- 
ments. It has also been questioned whether the second 
limit 77 — 0, physically associated to vanishing viscos- 
ity v — > 0, exists, or whether the limit is independent of 
the particular form of dissipation, such as hyperviscos- 
ity vs. ordinary viscosity. In a simple shell dynamical 
model, it was found by Leveque and She || that, while 
the limit v — v appeared to exist, the values of the scal- 
ing exponents in the putative inertial range depended 
upon the degree of the hyperviscosity. Subsequent nu- 
merical study of Navier-Stokes turbulence M found no 
similar dependence, although the low Reynolds number 
of the simulation makes such a result preliminary. It may 



be taken as the current "conventional wisdom" that the 
limit 77 — > exists and results in no ambiguity, while the 
limit L — > 00 does not exist (at least for moments) and is 
a source of anomalous scaling, if not an outright break- 
down of universality. However, the question which of the 
two limits exists, if either, and also the uniqueness of the 
limits must be regarded as open at the level of rigorous 
results. The issue frequently generates heated debates 
(e.g. 1). 

A simpler situation than the Navier-Stokes fluid for 
which many of the same questions arise is the problem of 
the turbulently convected passive scalar. The dynamics 
of the scalar field 9(r,t) is given by the linear equation 



(dt + v- V r )0 = k A r 9 + f 



(1) 



in which v(r,i) is the turbulent velocity and /(r, i) is a 
scalar source at macroscale L. The closest analogy exists 
for the so-called inertial- convective range of the passive 
scalar, in which the advecting velocities exhibit inertial- 
range scaling and the molecular diffusivity k is small, so 
that the scalar dynamics itself is dominated by convec- 
tion. In other words, this range corresponds to the ideal- 
ized limit of infinite Peclet number, Pc = UL/k = 00 (U 
is the typical velocity fluctuation at the macro length- 
scale). A dimensional theory for this situation analogous 
to that of Kolmogorov was developed soon after his by 
Obukhov and by Corrsin Q. However, similar evidence 
has been found to suggest that the Obukhov-Corrsin the- 
ory is flawed in the same way as K41 H|: while the exis- 
tence of the k —> limit is not threatened by the present 
data, the experiments suggest that there is a nontrivial 
dependence upon L due to intermittent cascade of the 
scalar from the macroscale. 

Recently, this "conventional picture" has obtained 
some support in a model of the passive scalar, first con- 
sidered by Kraichnan in 1968 @. In this model, the (in- 
compressible) velocity v and the source / are both taken 
to be Gaussian random fields delta-correlated in time. 
A regime mimicking the inertial-convective range of true 
turbulence is produced when the velocity covariance in 
space is taken to obey short-distance scaling: 



(« i (r)« i (0))~y *«-i3T C .[* < 

c 



(2) 



for r <C L. Kraichnan's "rapid-change" model is exactly 
soluble in the sense that there is no closure problem, i.e. 



1 



the Nth correlator of the scalar, 0jv(ri, ...,rjv;i) obeys 
a linear equation 



-H N Q N + F(r n ,r m )x 



pairs {nm} 



(3) 



which involves only itself and lower order correlators. 
The linear operator TLn is determined by the velocity 
statistics (see Eq.(^) below), while the inhomogeneous 
term involves the source covariance F. The new work 
has involved either (i) a physically-motivated ansatz for 
dissipative terms || or (ii) perturbative exploration of 
limiting regimes: C < 1 in || and \ < 1 in |Q. These 
studies all confirm the "conventional picture." Most re- 
cently |l2[ , the ^--expansion has predicted an interest- 
ing non-universality of the convective-range scaling ex- 
ponents, through dependence on the correlation time of 
the random velocity field (when that time is sufficiently 
small). In addition, the above theories all make detailed 
quantitative predictions of convective range scaling ex- 
ponents. Unfortunately, the two distinct approaches do 
not agree in detail and neither has complete control of er- 
rors. Numerical simulations |) 13 have not yet been able 
to probe the small £ or large d ranges where conflicting 
predictions occur. Thus, the results, while exciting, are 
still unconfirmed. 

We give a very direct, nonperturbative proof of the 
correctness of the "conventional picture" of inertial- 
convective range universality in the Kraichnan model. 
The full proofs are presented elsewhere [Q , but the most 
significant details may be easily explained here. The 
main step in our argument is to show that the station- 
ary solution of the model Eq.(|^) is unique for the ideal- 
ized k = case, when a certain regularity requirement 
is imposed. Our proof is inspired by an analogy of the 
model with the quantum many-body problem. In fact, 
the operator TLn for k = has the form of a quantum ki- 
netic energy operator in Nd dimensions with a position- 
dependent mass matrix M(R), or TLn — (2M(R)) _1 : 
PP. The mass matrix is defined in terms of the velocity 
covariance as 



[M(R)- 1 ] 



(vi{r n )vj(r m )}, 



(4) 



where R = (ri,...,rjy) is the ./V-particle configuration 
point and i,j ~ l,...,d, n,m — l,...,N. M(R) is there- 
fore a nonnegative matrix for every R. Because of the 
quantum-mechanical analogy it is very natural to look for 
solutions among square- integrable functions (Oat G £ 2 ). 
However, the perturbation studies [jllj suggest that the 
decay at large R is too slow for the solutions to be in C 2 . 
Since we are only interested anyway in the short-distance 
behavior, we confine the system to a bounded domain f2 
in d-dimensions, to remove this difficulty. On the bound- 
ary of Q we impose a fixed value To of the scalar, as would 
be appropriate in a real experiment with a temperature 



field T held fixed to To by thermal contact at the channel 
wall. Without loss of generality, we may take the bound- 
ary condition do = (Dirichlet b.c), by considering the 
fluctuation field 9 = T — T (as anticipated in Eq. ([!]).) 
We expect that our results will carry over also to more 
artificial geometries, e.g. the periodic b.c. used in nu- 
merical simulations |^| , although some parts of our proof 
below depend upon the Dirichlet b.c. 

There is a very simple formal argument which sug- 
gests the uniqueness of the k = solutions in C 2 . If we 
take Gn to be the inhomogeneous term in Eq.(|^), then 
the stationary solution of the equation should be given 
uniquely by 



©at — TL^Gn- 



(5) 



To permit this formal inversion, two things are required. 
First, it must be the case that TLn has no zero-modes. 
Secondly, it must also be shown that the range of TLn is 
closed (Fredholm condition). The insufficiency of the first 
condition alone may be explained more physically by ob- 
serving that the inversion is still ill-defined if the operator 
TLn has no zero-modes but if is a point of accumula- 
tion of the spectrum. In that case, the inverse TL N X is 
not a bounded operator and the righthand side of Eq. (||) 
need not exist. The Fredholm condition is guaranteed, in 
particular, if the "ground-state energy" of Hat is strictly 
positive. Unfortunately, there is a real danger that the 
condition is violated, because the mass matrix obviously 
has infinite mass eigenvalues at certain points R. This 
certainly happens whenever two points coincide, r„ = r m 
for n ^ m, since then the velocity covariance is degen- 
erate. Such large masses lead to the possibility of small, 
positive energies arbitrarily close to 0. Nevertheless, we 
have shown jl4j not only that zero-modes do not exist 
with Dirichlet b.c. but also that the "ground-state en- 
ergy" of Wat is then positive. We may note, incidentally, 
that a similar situation exists for periodic b.c. although 
(trivial) zero-modes certainly occur there, the constants. 
An equation may be developed for the cumulant part of 
the correlators, 6^, of the form TLn^n = in which 
the inhomogeneous term G C N is defined in terms of lower- 
order cumulants and is explicitly orthogonal to constants 
in the C 2 inner product. See Eq.(4.34) in |l5j| . One 
therefore has a similar formal solution <d N = H N 1 G C N . 
As above, this is indeed the unique solution if (i) con- 
stants are the only zero-modes of TLn and (ii) there is a 
positive spectral gap above the eigenvalue. We prove 
below that (i) holds, but (ii) is still open for periodic b.c. 

Establishing these facts requires a careful study of the 
mass matrix, carried out in (l4|] . It is shown there that 
zero eigenvalues of the inverse (Gramian) matrix G(R) = 
M(R) _1 occur only when some points r„,r m coincide, 
or are "fused." This result holds whenever the velocity 
spectrum in Fourier space is everywhere strictly positive. 
In general, the bad set of configurations R where zero 
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eigenvalues occur has K "clusters" of coinciding particles 
with Nk points in the fcth cluster, k = 1, K. It is shown 
in that the precise number of zero eigenvalues at 
each "bad" point is 52 fc=1 (iVfc — Furthermore, a lower 
bound is derived for the least eigenvalue v(R) of G(R) 
valid for all R, in terms of the minimum distance between 



any pair, p(R) = min. 



(Note that 



r m is the distance between pairs r„,r m . With periodic 



b.c. 



should be taken to be the least distance between 



all of their periodic images). The estimate that is proved, 
for < C < 2, is 



2/(R) > (const.) [^(R) 



(6) 



Establishing this bound is the crucial part of the proof. 
Its validity is suggested by a simple perturbation the- 
ory calculation in [p(R)]^ using the asymptotic short- 
distance formula (|^) for G(R). What is required is an 
argument that the first-order term in this expansion is 
non- vanishing. The difficulty is that there is a possible 
hierarchy of pair-separations {r nm } between the box size 
L and the minimum length p(R). However, a careful 
examination of cases using the general first-order per- 
turbation formula and the short-distance expression (j|) 
show that the leading term is always the first-order one. 
This argument yields Eq.(^). 

It is not hard to show from the above results for M(R) 
that constants are the only zero-modes of TIn- In fact, 
the stochastic representation holds that (@n, 'HnQn) = 
± fdR/dVP R (V)|(V.V R )eAr(R)| 2 , with P R (V) the 
multivariate Gaussian with covariance G(R) for the Nd- 
vector V. Because this probability density is strictly pos- 
itive for every R, V except for the "bad" points R of 
zero measure, it easily follows that any zero-mode satis- 
fies Vr0jv(R) = almost everywhere. Hence, all zero- 
modes are constants. In the case of Dirichlet b.c. this 
rules out existence of any zero-modes at all. 

The proof that the "ground-state energy" is strictly 
positive for Dirichlet b.c. uses the lower bound Eq.(||), 
giving 

(G N ,H N Q N ) > (const.) J dR [p(R)]^|V R e J v(R)| 2 (7) 

For any function g(R) whose Laplacian A R <? is positive 
and finite except for a singularity manifold T of codimen- 
sion > 2, the inequality holds that 

J dR | A R g(R)||<ys(R)| 2 

< 4 J dR\ A R .g(R)|- 1 |V R3 (R)| 2 |V R ^(R)| 2 , (8) 

with <p(R) any smooth function which = on the bound- 
ary dfl. This is the result of an integration by parts 
and the Cauchy- Schwartz inequality |l4|,[l(|. It may be 
applied to g(R) = [p(R)] c since A r5 (R) = 2((d - 



r y)[p(R)]~' y with 7 = 2 — £, which is positive as required 
when C > and d > 2. The Laplacian of g(R) has sin- 
gularities precisely on the "bad" set of points R with 
codimension d > 2 where mass eigenvalues become infi- 
nite. Thus, taking g(R) = [ ( o(R)] c and tp(R) = N (R) 
one obtains finally, if ( > 0, 

(d- -y) 2 

(Q N ,H N @ N ) > (const.) v - /; x 



J dR [p(R)me w (R)| 2 (9) 



for a set of Oat(R) dense in C 2 . Because p(R) < L 
everywhere, this implies that for 7 > 



(SjvjHtvQw) > (const.) 



■i-T|ie 



N 



(10) 



The last equation implies immediately a lower bound on 
the spectrum ~ (const. )L~ 7 whenever < ( < 2 and 
d > 2. In fact, by a criterion of Friedrichs, it can be 
shown from the two estimates Eqs.(^j) and (^) that the 
entire spectrum of Hn is discrete. See pq . Thus, the 
solution at k — is given uniquely in £r~hy Eq.(||) for 
Dirichlet b.c. Unfortunately, this proof does not work for 
periodic b.c. because the function g(R) has additional 
singularities on the codimension- 1 set of points where 
Tnm = r n r m , with r„» a periodic image of r„. These 
singularities contribute a surface term in the integration 
by parts argument which vitiates the derivation of the 
inequality Eq.(^). Despite the failure of this simple proof, 
we expect that there is still a gap in the spectrum above 
for periodic b.c. and all our results remain true. 

The proof of a unique solution in C 2 for n > can 
be given along the same lines but is much easier, since 
there is then obviously a spectral gap of size ~ ■ This 
solution formally corresponds to a state of zero Prandtl 
number, Pr = i> / k = 0. At high-wavenumbers it exhibits 
an infinitely long "incrtial-diffusive range" whose scalar 
spectrum decays with a large power. Such a range was 
proposed by Batchelor et al. in I|l7| and recently stud- 
ied in the "rapid-change" model |18| . This range follows 
after a lower wavenumber inertial-convective range a la 
Obukhov-Corrsin. Subsequently, one may consider the 
infinite Pcclct number limit, in which the high wavenum- 
ber end of the inertial-convective range will move off to 
plus infinity. It is not hard to show that limits of the 
positive-K solutions exist along subsequences for k — > 
and are C 2 solutions of the k — equations |1J]. Us- 
ing the uniqueness of the n — solutions, it then fol- 
lows that the k = solution is the limit of the positive- 
k solutions. It therefore corresponds to an idealized 
inertial-convective range of infinite extent. The same ar- 
gument applies to a broad class of dissipation terms be- 
sides the standard diffusion, including "hyperdiffusions" 
K p (—A r ) p 9. The conclusion is obtained that, for all these 
dissipation mechanisms, the limits exist in C 2 and coin- 
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cide with the unique n = solution (and therefore with 
each other). Universality is thus established. 

There is a simple physical interpretation of our re- 
sults which is worth emphasizing. Even in the limit as 
molecular diffusivity vanishes, n — > 0, there is an effec- 
tive diffusivity, or so-called eddy- diffusivity, acting on the 
mean scalar statistics which is generated through advec- 
tion by the random velocity. For iV-point statistics with 
N > 2 the effective diffusivity is scale-dependent, with 
Keddy(r) ~ D • at length-scale r. This was first ob- 
served for 2-point statistics (pairs of Lagrangian trac- 
ers) by Richardson p^| . He postulated a diffusion equa- 
tion of the same form as Eq.(p|) with an eddy-diffusivity 
tensor ify(r) scaling as Eq.(g) for £ = |. The matrix 
G(R) is the ./V-point generalization of the Richardson 
tensor and the Eq.@ for its lowest eigenvalue is just 
an expression of the expected scaling for K e ddy( r )- Note 
that the L-dependence of the spectral gap for k = is 
a direct consequence, since it should plausibly scale as 

^cddy(-fc) 



L -7 . It is important to observe that these 
eddy-diffusivity effects are only demonstrated for statis- 
tical correlations and not at the level of individual real- 
izations of the scalar field. This explains the paradoxical- 
looking fact that as k — ► the b.c. are satisfied by the 
statistical correlations, solving Eq.(|J), but not necessar- 
ily by the individual scalar fields, with dynamics Eq. ([[]). 
Even if the scalar source / is removed, it may not be 
possible to maintain the Dirichlet b.c. on 9(r, t) if k = 0. 
For such ideal dynamics the scalar field (when smooth) 
is simply transported unchanged along Lagrangian char- 
acteristics. Any point r in the interior of f2 with initial 
do (r) ^ which subsequently flows into the boundary dft 
will violate the b.c. at that point. For n > the b.c. on 
9 are maintained by the action of the molecular diffusion 
on the strong gradients of scalar concentration formed by 
such advcction to the boundary. This results in a thin 
"diffusive boundary layer" , with thickness vanishing as 
k — ► 0, across which the scalar intensity drops rapidly to 
zero. For our model, the b.c. for statistical correlations 
are maintained in the same way by the eddy-diffusivity, 
whose action creates a "turbulent boundary layer" across 
which iV-point functions drop to zero. 

The conclusions obtained here hopefully extend to the 
true passive scalar and, even more optimistically, to the 
Navier-Stokes fluid velocity. It is, at least, satisfying to 
have one model in which the "conventional wisdom" for 
these turbulent systems may be verified without question. 
Our analysis has some further interest in terms of the per- 
turbative studies fic||ll| . Those works obtained nontriv- 
ial zero-modes of local approximations to H.n by means 
of a Rayleigh-Schrodinger perturbation theory. This does 
not contradict our result that all zero-modes are con- 
stants, because we consider the global operator. Our re- 
sult just states that ambiguity in the perturbative treat- 
ments are removed by matching to the macro-scale L. 



In fact, our results confirm one of the basic assumptions 
used in the perturbative treatments: explicit matching to 
the dissipation scale may be replaced by a requirement 
of short-distance regularity. Even the weak requirement 
of local square-integrability is sufficient to select the so- 
lution properly matched to the dissipation range. 
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